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We study the three-dimensional structure formation when atoms are deposited onto a substrate 
with a decagonal quasicrystalline order. Molecular-dynamics calculations show that the adsorbate 
layer consists of ordered nano-scale domains with orientations determined by the underlying sub- 
strate symmetry. Depending on the relative strength of the interactions of adsorbate atoms with 
each other and with the substrate atoms, different morphologies are observed ranging from layer- 
by-layer growth to cluster formation. We also find that the film thickness likewise affects the overall 
structure of the growing film: Depending on the relative strength of the interaction between ad- 
sorbate atoms, a structural transition of the configuration of the adsorbate layers closest to the 
substrate can occur as the number of deposited layers increases. 



Epitaxial growth produces structures with tailored 
properties, which do not exist in three dimensions (3D) 
due to thermodynamical restrictions. For example, it 
allows for the formation of nano-sized structures and 
thereby opens up new perspectives in science and tech- 
nology. The large amount of experimental results that 
have been obtained on self-size selecting growth over the 
last years poses an immense challenge for theorists. 

In this context, quasicrystalline media are of particular 
theoretical and experimental interest, because they are 
neither disordered like an amorphous solid nor periodic 
like a crystal. Contrary to common vitreous structures, 
quasicrystals do show perfect long-range order. Their 
surfaces can routinely be prepared, they are increasingly 
used as substrates for hetero-epitaxial growth and have 
been studied intensively over the last yearsi^^i^i^i^. How- 
ever, theoretical studies of this kind are complicated by 
the fact that the precise 3D structures of quasicrystals, 
i.e., their atomic configurations, are not readily available. 

Fliickiger et al.^ have investigated the interface formed 
by an Al film grown on the surface of a decagonal 
AlyoCoisNiis (Al-Co-Ni) substrate and found that the 
mismatch at the interface causes the formation of well- 
oriented nanometric single crystals of Al, satisfying the 
epitaxial conditions on a local scale. Total energy cal- 
culations based on the simplifying assumption that the 
deposited layer interacts as a rigid lattice with the un- 
derlying substrate satisfactorily reproduced their exper- 
imental findings. In subsequent experimental work by 
Ferralis et al.-, the adsorption of Xe on the surface of 
decagonal Al-Co-Ni was investigated and results quali- 
tatively similar to those found by Fliickiger et alfi were 
obtained. 

Recently, physisorption of Xe on decagonal Al-Co-Ni 
has been simulated using 3D Grand Canonical Monte 
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Carlo techniques^i^. Setyawan et al^ have shown that 
single locally hexagonal domains align only along fixed 
directions that have fivefold symmetry. However, since 
the size of the unit cell used is of the order of the domain 
size, it is difficult to observe the formation of multiple 
domains. This circumstance points towards the need 
of using much larger substrate sizes in the simulations 
that, additionally, will also reduce the effect of the sam- 
ple boundaries. 

For the systems we consider, i.e., adsorption of atoms 
with an essentially isotropic potential on a quasicrys- 
talline substrate and not showing effects like intermix- 
ing, real-space experimental data on the formation of 
domains have been lacking. The dependence of the do- 
main size and shape on experimental parameters is like- 
wise not known. In a recent 2D study- we have pointed 
towards the importance of the relative strength of the 
mutual interaction of adsorbate atoms (adatoms) com- 
pared to their interactions with the substrate atoms, 
which we quantified by the parameter 77. We found that 
with increasing 77 the adatoms assemble into well-formed 
and orientationally locked crystalline domains, qualita- 
tively similar to the experimental findings^. However, 
the growth of adsorbate layers (adlayers) on a quasicrys- 
tal is a 3D process and should be treated as such. 

In this article, we study the adlayer growth by carry- 
ing out extensive numerical simulations of the 3D depo- 
sition of 3.5 monolayers (ML) of 11 000 adatoms. As a 
substrate, we used 6 000 atoms in a surface bilayer of di- 
mensions 150 X 150 X 2.05 A3 which was extracted from a 
model structure of the decagonal Al-Co-Ni quasicrystal 
recently proposed by Deloudi and Steurer— . The large 
lateral size of the unit cell permits us to observe the for- 
mation of multiple domains and reduces artefacts gener- 
ated by boundary conditions. We find that in addition to 
77, the film thickness also plays a role in determining the 
overall structure of the growing film, as one expects. The 
structure obtained after the formation of the first layer 
therefore does not necessarily determine the overall struc- 
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ture of the final film: As more and more layers are added 
to the film, the influence of the substrate on the first layer 
can diminish in favor of its natural structure. Thus, with 
increasing number of layers the first layer can undergo a 
structural transition. We also find that the parameter 77 
strongly influences the mode of growth: Layer-by-layer 
growth at ?7 = f gradually crossing over into 3D growth 
(clustering) at around ry = 2. For rj > 1.35 we observe 
that the adatoms self-assemble into locally well-ordered 
domains with orientations determined by the symmetry 
of the underlying substrate, similar to what is observed 
in the experiments^'^, and confirming our earlier 2D nu- 
merical modeling and simulations^. 

The adatoms are deposited by injecting them at an 
initial height zq above the substrate surface with a uni- 
form initial velocity vq. We assume that the substrate 
remains rigid and use a cooling mechanism to maintain 
the adlayer at a constant temperture T, as we will ex- 
plain below. The interaction between adatoms as well as 
between adatoms and substrate atoms is assumed to be 
pairwise and of Lennard-Jones type, with the character- 
istic energy and length scales given by e and a, respec- 
tively. 

Letting and ri denote the coordinates of a substrate 
atom a and an adatom i, respectively, we numerically 
evaluate the equations of motion 
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where the parameter rj is used to adjust the relative 
strength of the adatom interactions with respect to the 
substrate. 

The adatoms are maintained at a constant temper- 
ature T by implementing a thermostat proposed by 
Gilmore and Sprague^*^: Excess kinetic energy of the 
adatoms is dissipated by an additional friction force (T) 
that is given by the product of a damping constant 7 
times the excess kinetic energy with respect to the de- 
sired average kinetic energy. 

The numerical simulations are carried out in a unit cell 
of 150 X 150 X 50 A"^ with periodic boundary conditions in 
the lateral {x, y) directions and a hard ceiling at 2; = 50 A. 
The quasicrystalline substrate is taken as a bilayer with 
a layer spacing of 2.05 A containing about 3 000 atoms in 
each layer. It is assumed that the substrate is sufficiently 
large and the effects of artificial periodicity induced by 
the boundary conditions can be neglected. For the inter- 
action parameters we used a — 2.55 A corresponding to 
the lattice parameter of Al metal, while e is set to 0.25 
eV for practical purposes. 

Adatoms were inserted in batches of 180 from a height 
of zo = 25 with an initial velocity of wo = 5 (in time units 
where the mass of the adatoms is set to unity) . In order 
to avoid correlations between the incident adatoms they 
were injected at random lateral {x, y) positions, such that 
the distance between any two of them was at least 3fT. 
The substrate was maintained at a temperature kT = 0.1 



eV, and a damping coefficient 7 = 0.25 was used. It was 
determined in advance that a waiting time of ^wait = 25 
time units was sufficient for the incident adatoms to dissi- 
pate their energy to the adlayer, and for the latter in turn 
to relax back to the thermostat temperature T. Thus, in 
our simulations we inject 180 adatoms at a time, wait 
for twait before inserting another batch. The equations 
of motion, Eqs. ([T]), were integrated using a fourth-order 
Runge-Kutta method with adaptive step sizeii. 

We first investigate the occupancy of the adlayers as a 
function of the total number of adatoms deposited onto 
the substrate. For 77 values of 1, 1.5, and 2, adlayer struc- 
tures of 3.5 ML containing 10980 adatoms were grown, as 
described above. For each 77 the individual layer heights 
were determined from the height distributions of the 
adatoms. The results are shown in Fig. [2 where each 
continuous curve consists of 61 data points, correspond- 
ing to the 61 batches of 180 adatoms injected upon the 
surface. We see that for ry values close to one the de- 
position curves are rather steep with a slope of nearly 
one, reflecting a layer-by-layer type of growth. For larger 
•q values the characteristic deposition curves are flatter, 
implying that multiple layers are grown simultaneously. 

In order to investigate the influence of the substrate 
structure on the growing film as a function of film thick- 
ness, we follow the structure of the adlayer closest to 
the substrate. Fig. [2] displays a set of intensity patterns 
obtained from the Fourier transformation of the adatom 
structure of the first layer for -q = 1 (left column) and 
1] = 2 (right column) for the deposition of 1.5, 2.5, and 
3.5 ML of adatoms (from top to bottom). Note that for 
r] = 1 the interface layer of the growing film shows a 
tenfold symmetry, induced by the quasicrystalline sub- 
strate, that does not change appreciably with increasing 
film thickness. This is in contrast to the growth mode 
at 77 = 2, where the structure of the film at the interface 
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FIG. 1: (color online) Layer occupancies in layers L = 1 — 4 as 
a function of the total number of adatoms deposited. Shown 
are results for relative interaction strengths rj = 1, 1.5, and 2. 
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FIG. 2: Intensity patterns obtained from the Fourier trans- 
formation of the adatom structure of the first layer for — 1 
(left column) and rj = 2 (right column) for the deposition of 
1.5, 2.5, and 3.5 ML of adatoms (from top to bottom). In (f), 
one set of six spots is marked by a hexagon to emphasize the 
sixfold symmetry of the domain that gives rise to the pattern. 
The azimuthal orientation of the next domain is 12° away. 




FIG. 3: (color online) Overall adatom configuration at 3.5 ML 
coverage for = 1 (a) and 77 = 2 (b). Atoms belonging to 
a particular layer have the same color. The interface layer 
{L — 1) is yellow. 



changes from quasicrystalline to local hexagonal domains 
with orientations locked to the substrate symmetry. One 
sixfold-intensity pattern is highlighted by the hexagon in 
Fig. [2] (f). The layer thicknesses have been indicated in 
Fig. [T] as thin vertical lines. The layer occupancies corre- 
sponding to the patterns shown in Figs. [2] (a) - (f) have 
also been marked in Fig. [TJ 

We now turn to the real-space configuration of the film. 
Figs. [3] (a) and (b) present the overall atomic configura- 
tion of 3.5 ML film thickness for 77 = 1 and rj — 2, respec- 
tively. It is evident that rj — 1 results in a layer-by-layer 
growth with relatively high occupancies of the deeper lay- 
ers, while the film produced with ry = 2 is highly corru- 
gated. These results are consistent with the deposition 



curves presented in Fig. [T] The ry = 2 film displays well- 
formed hexagonal domain structure. In a previous work, 
we have used Voronoi tessalation to reveal the domains of 
hexagonal structure in a 2D configuration of adatoms^. 
We apply this technique to the configuration of adatoms 
in the interface layer calculated with rj — 2, Fig. [3] (b). 
The result is shown in Fig. 3] Voronoi polygons other 
than hexagons are colored dark (light violet for n > 7- 
gons and dark violet for n < 5-gons) , while hexagons have 
been colored according to their orientation relative to the 
X-axis as indicated by the color-coded scale from — 60° 
on the left-hand side of the figure. We distinguish at least 
four well-oriented hexagonal domains of orientations 6° 
(dark blue), 18° (light blue), 42° (light green), and 54° 
(orange), consistent with the corresponding Fourier in- 
tensity pattern in Fig. [5] (f). The domains have sharp 
boundaries consisting of non-hexagonal polygons, i.e., 
topological defects. We also see some continuous transi- 
tions between domains with orientations differing by 12°. 
Notice that for larger differences of orientations of neigh- 
boring domains the domain boundaries seem to consist 
mostly of topological defects. Also note the large do- 
main sizes of the order of 50 A in Fig. [H This shows di- 
rectly the necessity of using appreciably large substrates 
in order to accomodate the formation and interactions of 
multiple domains while reducing the effects of boundary 
conditions. 
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FIG. 4: (color online) Voronoi tessalation of the real-space 
configuration of adatoms (black spots) in the interface layer. 
Hexagons have been color coded with respect to their orien- 
tation relative to the a;-axis. Voronoi n-gons with n > 7 and 
n < 5 have been colored in light and dark violet, respectively, 
in order to clearly mark non-hexagonal polygons. 

Our 3D results on a hetero-epitaxial system consiting 
of a quasicrystal and an ordinary crystal show that the 
structural misfit can lead to the formation of nano-scale 
crystalline domains. The aperiodic substrate further acts 
as a structural template to azimuthally align these do- 
mains along the distinct symmetry directions of the qua- 
sicrystal. Moreover, we find that the domain formation 
and structure within a given adatom layer is consistent 
with our earlier findings on a 2D system in which single 
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layers of adatoms were confined to a plane^ . These results 
are also in agreement with experimental findings2,!^i^. 

The clustering growth resulting in a well-formed local 
crystalline structure at 77 = 2, as clearly seen in Fig.[3](b), 
shows that the formation of the growing film is dominated 
by the relatively strong adatom-adatom interactions in 
comparison with the adatom-substrate interaction which 
causes the injected adatoms to preferentially coagulate 
with already present and cooled adatoms forming nano- 
clusters that serve as seeds. This is in contrast to the 
wetting behavior of the growing film at 77 = 1, shown 
in Fig. [3] (a). Hence, by varying the parameter rj it is 
possible to continuously cover a wide range of growth 
modes. For ?; = 2, we also observe a structural transition 
in the interface layer of the film as it gradually grows. 
The presence of multiple adlayers on top of the interface 
layer shows a strong tendency to revert the latter to its 
natural sixfold orientational symmetry. We believe that 
this transition is caused by the structural stiffening of the 
film with increasing number of adlayers. 

Experimental results of Fliickiger et al.^ for Al films 
grown on Al-Co-Ni show a transition from quasicrys- 
talline structure at 1 ML to a structure of well-oriented 
nanoscale single crystals of Al at coverages of 2 ML and 
above (see also Ref.— for similar results of Xe adsorbed 
on Al-Co-Ni). We observe qualitatively similar results for 
•q = 2. In fact, additional simulations not presented here 



show that such a transition persists for values of 77 as low 
as 1.35. It is very likely that such a threshold value de- 
pends also on the damping coefficient which determines 
the rate at which arriving adatoms cool down on the 
adlayer. A lower damping constant will increase the re- 
laxation times for these atoms and could thus permit the 
exploration of a larger set of energetically more favorable 
local configurations. Although the layer-by-layer growth 
seen for 77 ~ 1 suggests a less pronounced dependency in 
this regime, it is possible that the structure of the wetting 
adlayer might nevertheless be affected by larger thermal 
relaxation times. We observed in a previous simulated 
annealing study of a single adlayer a dependence of the 
local structure on cooling rate^. 

The present article considers the formation of struc- 
ture on a quasicrystalline substrate in its dependence 
on the relative strength of adatom interactions and film 
thickness. A study investigating the effect of substrate 
temperature as well as the thermal relaxation rate is cur- 
rently in progress. 
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